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We report recent results of the Dalitz plot analyses of D and Ds decays performed by the BABAR 
collaboration, and point out some of the important applications of these results. 
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INTRODUCTION 

The amplitudes describing D and Ds meson weak de- 
cays into final states with three pseudo-sealers are dom- 
inated by intermediate resonances that lead to highly 
nonuniform intensity distributions in the available phase 
space. The results of the Dalitz plot analysis of these 
decays are playing increasingly important role in flavor 
physics, particularly in the extraction of the CP- violating 
phase 7 = arg(-14^F4/y^^V;;) of the quark mixing 
{i.e., CKM) matrix by exploiting interference structure 
in the D Dalitz plot from the decay B^ -^ DK^ |l|] and 
in the measurement of D^-D^ mixing parameters. 



DETECTOR 

We perform these analyses using e^e~ collision data 
collected at and around 10.58 GeV center-of-mass (CM) 
energy with the EaBAR detector [1] at the PEP-II stor- 
age ring. Tracking of charged particles is provided by 
silicon detector and a drift chamber operating in a 1.5- 
T magnetic field. Particle types are identified using 
specific ionization energy loss measurements in the two 
tracking devices and Cherenkov photons detected in a 
ring-imaging detector. The energy of photons and elec- 
trons is measured with an electromagnetic calorimeter. 
In case of neutral £>-meson decays, we distinguish D^ 
from D^ by reconstructing the decays £>*+ — > D°7r+ and 
D*^ -^ D^Ti^ . For each decay mode, we estimate the 
signal efficiency as a function of position in the Dalitz 
plot using simulated signal events generated uniformly in 
the available phase space, subjected to the same recon- 
struction procedure applied to the data, and corrected 
for differences in particle-identification rates in data and 
simulation. 



final state can be characterized as a coherent sum of all 
relevant quasi-two-body D/Ds ^ (r ^ AB)C isobar 
model resonances, A = Ylr^re^'^'' Ar{s). Here s = m\g, 
and Ar is the resonance amplitude. We obtain the coef- 
ficients Or and 4>r from a likelihood fit. The probability 
density function for signal events is |^|". 

Unless stated otherwise, for S-, P-, and D-wave (spin 
= 0, 1, and 2, respectively) resonant states we use the 
Breit-Wigner amplitude: 
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where Mq (Fq) is the resonance mass (width) [3|, L is 
the angular momentum quantum number, p is the mo- 
mentum of either daughter in the resonance rest frame, 
and po is the value of p when s = Mq. The func- 
tion J^L is the Blatt-Weisskopf barrier factor [J]: .Fq = 
and T2 ~ 1/ 



1, Ti = 1/V1 + Rp, and T2 = l/V9 + 3Rp^ + Rp^, 
where we take the meson radial parameter R to be 
1.5 GcV~^. The quantity Ml is the spin part of 
the amplitude: A^o = constant, A^i ex ^2pA.pci ^^^id 
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where pi is the 3- 
momentum of particle i in the resonance rest frame. 
The fit fraction for a resonant process r is defined as 
fr = f \arAr\'^ dr/ J \A\'^ dr, where dr is a phase-space ele- 
ment. Due to interference among the contributing ampli- 
tudes, the fr do not sum to one in general. In all cases, 
we model small incoherent background empirically from 
data. 



ANGULAR MOMENTS 



DALITZ PLOT PARAMETRIZATION 

The complex quantum mechanical amplitude A that 
describes decays to three particles A, B and C in the 
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For D and Ds decays to three spinlcss particles, the 
Dalitz plot uniquely represents the kinematics of the final 
state. The angular distributions provide further informa- 
tion on the detailed event-density variations in various 
regions of the phase space in a different form. We define 
the helicity angle Oh for decays D° ^ (r ^ AB)C as the 
angle between the momentum of A in the AB rest frame 
and the momentum of AB in D^ rest frame. The mo- 
ments of the cosine of the helicity angle, Y^* {cos Oh), are 



defined as the efficiency-corrected invariant mass distri- 
butions of events when weighted by spherical harmonic 
functions 
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where m is the invariant mass of the AB system and the 
Pi are Legendre polynomials of order /: 



Pl{x) Pn{x) dx = Sin- 
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These angular moments have an obvious physical sig- 
nificance. Since spherical harmonic functions are the 
eigen- functions of the angular momentum, the Dalitz plot 
of a three-body decay can be represented by the sum of 
an infinite number of spherical harmonic moments in any 
two-body channel. In a region of the Dalitz plot where 
S- and P- waves in a single channel dominate, their am- 
plitudes are given by the following Legendre polynomial 
moments, 
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where |5| and \P\ are, respectively, the magnitudes of 
the S- and f-wave amplitudes, and 9sp = Os — (^p is the 
relative phase between them. It is worth noting that this 
partial-wave analysis is valid, in the absence of higher 
spin states, only if no interference occurs from the cross- 
ing channels. 



DALITZ PLOT ANALYSIS OF D° -^ A A' 
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FIG. 1: Dalitz plot for D° -^ A'A+Tr" [3 data (a), and the 
corresponding squared invariant mass projections (b-d). In 
plots (b-d), the dots with error bars are data points and the 
solid lines correspond to the best isobar fit models. 
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FIG. 2: LASS (solid line) and E-791 (dots with error bars) 
Ktt S-wave amplitude (a) and phase (b). The double headed 
arrow indicates the mass range available in D^' — > A^A'*'7r°. 



The A:±7r° systems from the decay D" -^ A'" A'+7r° ^ 
can provide information on the Kit S-wave amplitude in 
the mass range 0.6-1.4 GeV/c^, and hence on the pos- 
sible existence of the k(800), reported to date only in 
the neutral state {k'^ -^ K^n^) [2|. If the k has isospin 
1/2, it should be observable also in the charged states. 
Results of the present analysis can be an input for ex- 
tracting the CKM phase 7 by exploiting interference in 
the Dalitz plot from the decay B± -^ P>](-k+-k°^^ 0- 

We perform the analysis on 385 fb~^ data using the 
same event-selection criteria as in our measurement of 
the branching ratio of the decay D° -^ K'K+n^ 0. To 
minimize uncertainty from background shape, we choose 
a high purity (^98%) sample using 1855 < nioa < 
1875 MeV/c^ and find 11278 ± 110 signal events. The 
Dalitz plot for these events is shown in Fig. [IJa). 
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FIG. 3: The phase-space-corrected A^ A^ S- and P-wave am- 
plitudes, |S| and \P\, respectively, (a) Lineshapes for (solid 
line, blue) /o(980), and (broken line, blue) ao(980). (b) Line- 
shape for (j){1020) (solid line, blue). In each plot, solid circles 
with error bars correspond to values obtained from the model- 
independent analysis. In (a), the open triangles (red) corre- 
spond to values obtained from the decay Z)" -^ K~ K^ K^. 
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FIG. 4; Legendre polynomials moments for the K^-k^ (columns I, II) and K K^ (columns III, IV) channels of £>" 
The circles with error bars are data points and the curves (red) are derived from the fit functions. 
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For 13° decays to A'^tt" S'-wavc states, we consider 
three amplitude models: LASS amplitude for K~'k'^ — > 
K~t:'^ elastic scattering [8|, |9|, the E-791 results for the 
K^TT'^ S'-wavc amplitude from a partial-wave analysis of 
the decay D+ -^ Ar~7r+7r+ [lO[, and a coherent sum of 
a uniform nonresonant term plus Breit-Wigner terms for 
k(800) and ArQ(1430) resonances. 

In Fig. [2] we compare the Kir S-wave amplitude from 
the E-791 analysis to the LASS amplitude. The 
LASS Ktt S-wave amplitude gives the best agreement 
with data and we use it in our nominal fits {x^ proba- 
bility 62%). The Ktt 5- wave modeled by the combina- 
tion of k(800) (with parameters taken from Ref. [6|]), a 
nonresonant term and i^Q(1430) has a smaller fit prob- 
ability (x^ probability < 5%). The best fit with this 
model (x^ probability 13%) yields a charged k of mass 
(870 ± 30) MeV/c2, and width (150 ± 20) MeV/c^, sig- 
nificantly different from those reported in Ref. [6[ for the 
neutral state. This does not support the hypothesis that 
production of a charged, scalar k is being observed. The 
E-791 amplitude [10[ describes the data well, except near 
threshold. We use it to estimate systematic uncertainty 
in our results. 

We describe the 13° decay to a K~ K^ 5- wave state 
by a coupled-channel Breit-Wigner amplitude for the 
/o(980) and ao(980) resonances, with their respective 
couplings to TTTT, KK and rjir, KK final states [9[. Only 
the high mass tails of /o(980) and ao(980) are observable, 
as shown in Fig. [3l 

We find that two different isobar models describe the 
data well. Both yield almost identical behavior in invari- 
ant mass (Fig. [T|d-[TJ1) and angular distribution (Fig. [4]). 
The dominance of D° -^ K*+K- over D° -^ K*- K+ 
suggests that, in tree-level diagrams, the form factor for 
D^ coupling to K*~ is suppressed compared to the corre- 
sponding K~ coupling. While the measured fit fraction 
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FIG. 5: Results of the partial- wave analysis of the K~ K^ 
system, (a) Cosine of relative phase Qsp = 6s — Op, (b) two so- 
lutions for 6sp, (c) P-wave phase for 0(1020), and (d) S'-wavc 
phase derived from the upper solution in (b). Solid bullets 
are data points, and open circles (blue) and open triangles 
(red) correspond, respectively, to isobar models I, II. 



for D^ -^ K*^K agrees well with a phenomenological 



prediction [ll| based on a large SU(3) symmetry break- 
ing, the corresponding results for Z)° — > K*~K'^ and the 
color-suppressed D^ -^ (J)tt° decays differ significantly. It 
appears from Table [J that the A'+tt*^ S-wave amplitude 
can absorb any AT* (1410) and /2(1525) if those are not 
in the model. The other components are quite well es- 
tablished, independent of the model. From Table HI the 
strong phase difference, Sd, between the D^ and D° de- 
cays to AT* (892)+ AT" state and their amplitude ratio, r^. 



TABLE I: The results obtained from the D" — > K K^tv^' DaUtz plot fit [3|. The errors are statistical and systematic, 
respectively. We show the ao(980) contribution, when it is included in place of the /o(980), in square brackets. 





Model I 


Model If 1 


State 


Amplitude, a^ 


Phase, (f>r (°) 


Fraction, fr (%) 


Amplitude, a^ 


Phase, (t>,. (°) 


Fraction, fr (%) 


7^* (892)+ 


1.0 (fixed) 


0.0 (fixed) 


45.2±0.8±0.6 


1.0 (fixed) 


0.0 (fixed) 


44.4±0.8±0.6 


A'* (1410)+ 


2.29±0.37±0.20 


86.7±12.0±9.6 


3.7±1.1±1.1 








K+n^S) 


1.76±0.36±0.18 


-179.8±21.3±12.3 


16.3±3.4±2.1 


3.66±0.11±0.09 


-148.0±2.0±2.8 


71.1±3.7±1.9 


0(1020) 


0.69±0.01±0.02 


-20.7±13.6±9.3 


19.3±0.6±0.4 


0.70±0.01±0.02 


18.0±3.7±3.6 


19.4±0.6±0.5 


/o(980) 


0.51±0.07±0.04 


-177.5±13.7±8.6 


6.7±1.4±1.2 


0.64±0.04±0.03 


-60.8±2.5±3.0 


10.5±1.1±1.2 


[ao(980)°] 


[0.48±0.08±0.04] 


[-154.0±14.1±8.6] 


[6.0±1.8±1.2] 


[0.68±0.06±0.03] 


[-38.5±4.3±3.0] 


[11.0±1.5±1.2] 


/2(1525) 


1.11±0.38±0.28 


-18.7±19.3±13.6 


0.08±0.04±0.05 








7^*(892)- 


0.601±0.011±0.011 


-37.0±1.9±2.2 


16.0±0.8±0.6 


0.597±0.013±0.009 


-34.1±1.9±2.2 


15.9±0.7±0.6 


A'*(1410)- 


2.63±0.51±0.47 


-172.0±6.6±6.2 


4.8±1.8±1.2 








K-n^iS) 


0.70±0.27±0.24 


133.2±22.5±25.2 


2.7±1.4±0.8 


0.85±0.09±0.11 


108.4±7.8±8.9 


3.9±0.9±1.0 



are given by: So = -35.5° ± 1.9° (stat) ±2.2° (syst) and 
ro = 0.599 ± 0.013 (stat) ± 0.011 (syst) 0. Systematic 
uncertainties in quantities in Table [J arise from exper- 
imental effects {e.g., efficiency parameters, background 
shape, particle-identification), and also from uncertainty 
in the nature of the models used to describe the data 
(e.g., A'tt S'-wave amplitude and resonance parameters). 
We show the Legendre polynomials moments in Fig. 3] 
for the i^+7r° and K~K~^ channels, for ^ = 0— 7. We use 
the relations of Eq. [5] to evaluate \S\ and \P\ shown in 
Fig. [3l and 9sp shown in Fig. [5l for the K~K^ channel 
in the mass range mx-K+ < 1-15 GcV/c^. The mea- 
sured values of jS"! agree well with those obtained in the 
analysis of the decay D^ — s- K~K'^K^ [l^] and also with 
either the /o(980) or the ao(980) lineshape. The mea- 
sured values of |P| are consistent with a Brcit-Wigner 
lineshape for (f>{lQ2Q). 



DALITZ PLOT ANALYSIS OF D° ^ tt 7r+7r° 

An important component of the program to study 
CP violation is the measurement of the angle 7 of 
the unitarity triangle related to the Cabibbo-Kobayashi- 
Maskawa quark mixing matrix. The decays B -^ 
£)i*)0]^i*) Qg^jj j-,g used to measure 7 with essentially no 
hadronic uncertainties, exploiting interference between 
b — > ucs and b — > cus decay amplitudes. The most 
effective method to measure 7 has turned out to be 
the analysis of the D-decay Dalitz plot distribution in 
B^ -^ DK^ with multi-body D decays [H. This 
method has only been used with the Cabibbo-favored 
decay D ^ A'^+tt" 0, [H. We perform the first 



CP-violation study of B^ -^ DK^ using a multibody, 
Cabibbo-suppressed D decay, D -^ tt+tt^tt". 

We determine the parameters a^, (f^r, and fr by fit- 
ting a large sample of I?° and D^ mesons, flavor-tagged 
through their production in the decay D*+ -^ I?°7r+ [7|. 
Of the D candidates in the signal region 1848 < tti^jo < 
1880 McV/c^, we obtain from the fit 44780 ± 250 signal 
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FIG. 6: Dalitz plot and invariant mass-squared projections 
for the D" -^ tt^tt+tt^ decay excluding D" -> A'>°. 



and 830 ± 70 background events. 

Table HIl summarizes the results of this fit, with system- 
atic errors obtained by varying the masses and widths of 
the p(1700) and a resonances and the form factors, and 
also varying the signal efficiency parameters to account 
for uncertainties in reconstruction and particle identifi- 
cation. The Dalitz plot distribution of the data is shown 
in Fig. [6ja-d). The distribution is marked by three de- 
structively interfering pir amplitudes, suggesting a final 
state dominated by / = [l6|. We show the Legendre 
polynomials moments in Fig. [7] for the tt+tt^ and 7r^7r+ 
channels, for 1 = 0—7. The agreement between data 
and fit is again excellent. Unlike in case of the decay 
D^ -^ K^K'^tt'^, we cannot use the relations of Eq. [5]to 
evaluate \S\ and \P\, and 9sp in any of the two-body tttt 
channels because of the contributions from cross-channels 
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FIG. 7: Legendre polynomials moments for the tv'^tt^'' (columns I, II) and it tv^ (columns III, IV) channels of Z?" — » tv tt^tt". 
The circles with error bars are data points and the curves (red) are derived from the fit functions. 



in the entire available mass-range. 
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FIG. 8: (a) The D° -^ K'^tt'tv+ Dalitz distribution from 



D* 



D TT events, and projections on (b) m^ = 
- = "^KgTT-' a^nd (d) m^+^_. D° -^ Kl-K+Tv- from 

The curves are the 



D*^ —^ D TT^ events are also included 
model fit projections. 



DALITZ PLOT ANALYSIS OF D° -^ A'Stt+tt- 



The Dalitz plot analysis of the decay D° — > K'^tt^tt is 
also motivated by its application to the measurement of 
CKM phase 7 [i3|. We determine the £>° -^ K^tt+tt- de- 
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FIG. 9: The invariant mass distribution of the reconstructed 
Ds candidate in the decay D^ -^ K^ K~ -k^ . For the Dalitz 
plot analysis we use events in the mass window shown by 
vertical arrows. The results are preliminary. 



cay amplitude from an imbinned maximum-likelihood fit 
to the Dalitz plot distribution of a high-purity (~ 98%) 
D^ sample from 390328 D*^ -^ D^ii^ decays recon- 
structed in 270 fb^ of data, shown in Fig. [8] 

The decay amplitude is expressed as a coherent sum 
of two-body resonant terms and a uniform non-resonant 
contribution. For r = p(770) andj9(1450) we use the 
functional form suggested in Ref. [18| . while the remain- 
ing resonances are parameterized by a spin-dependent 
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FIG. 10; (a) The _D+ ^ iC+iC"7r+ Dalitz distribution, and 
projections on (b) m^+^_, (c) m^-_^+, and (d) m^+^+. 
The curves are the model fit projections. The results are 
preliminary. 



TABLE II; The results obtained from the D° ^ 7r~7r+7r° 
Dalitz plot fit [l[. The errors are statistical and systematic, 
respectively. We take the mass (width) of the a meson to be 
400 (600) MeV/c^ 



State 


a. (%) 


4>rn 


/.(%) 


p(770)+ 


100 





67.8±0.0±0.6 


p(770)° 


58.8±0.6±0.2 


16.2±0.6±0.4 


26.2±0.5±1.1 


p(770)- 


71.4±0.8±0.3 


-2.0±0.6±0.6 


34.6±0.8±0.3 


p(1450)+ 


21±6±13 


-146±18±24 


0.11±0.07±0.12 


p(1450)° 


33±6±4 


10±8±13 


0.30±0.11±0.07 


p(1450)- 


82±5±4 


16±3±3 


1.79±0.22±0.12 


p(1700)+ 


225±18±14 


-17±2±3 


4.1±0.7±0.7 


p(1700)° 


251±15±13 


-17±2±2 


5.0±0.6±1.0 


p(1700)- 


200±11±7 


-50±3±3 


3.2±0.4±0.6 


/o(980) 


1.50±0.12±0.17 


-59±5±4 


0.25±0.04±0.04 


/o(1370) 


6.3±0.9±0.9 


156±9±6 


0.37±0.11±0.09 


/o(1500) 


5.8±0.6±0.6 


12±9±4 


0.39±0.08±0.07 


/o(1710) 


11.2±1.4±1.7 


51±8±7 


0.31±0.07±0.08 


/2(1270) 


104±3±21 


-171±3±4 


1.32±0.08±0.10 


o-(400) 


6.9±0.6±1.2 


8±4±8 


0.82±0.10±0.10 


Non-Res 


57±7±8 


-11±4±2 


0.84±0.21±0.12 



relativistic Breit-Wigner distribution. The model con- 
sists of 13 resonances leading to 16 two-body decay am- 
plitudes and phases (sec Tabic ITTT|) . plus the non-resonant 
contribution, and accounts for efficiency variations across 
the Dalitz plane and the small background contribution. 
All the resonances considered in this model are well es- 
tablished except for the two scalar tttt resonances, u and 
cr', whose masses and widths are obtained from our sam- 
ple [l9[. Their addition to the model is motivated by an 
improvement in the description of the data. 



The possible absence of the a and cr' resonances is con- 
sidered in the evaluation of the systematic errors. In this 
respect, the K-matrix formalism [20| provides a direct 
way of imposing the unitarity constraint that is not guar- 
anteed in the case of the Breit-Wigner parametrization 
and is suited to the study of broad and overlapping res- 
onances in multi-channel decays. We use the K-matrix 
method to parameterize the tttt S-wave states, avoiding 
the need to introduce the two a scalars. A description of 
this alternative parametrization can be found in Ref . [21 1 . 



Component 


Re{are'f-} 


/m{a,e'*"-} 


/. (%) 


/s:*(892)- 


-1.223 ±0.011 


1.3461 ± 0.0096 


58.1 


A-o*(1430)- 


-1.698 ±0.022 


-0.576 ± 0.024 


6.7 


A-2*(1430)- 


-0.834 ±0.021 


0.931 ±0.022 


6.3 


A'*(1410)" 


-0.248 ± 0.038 


-0.108 ±0.031 


0.1 


A'*(1680)~ 


-1.285 ±0.014 


0.205 ±0.013 


0.6 


A-* (892)+ 


0.0997 ± 0.0036 


-0.1271 ±0.0034 


0.5 


7^0* (1430)+ 


-0.027 ±0.016 


-0.076 ±0.017 


0.0 


is:2*(1430)+ 


0.019 ±0.017 


0.177 ±0.018 


0.1 


p(770) 


1 





21.6 


Lj(782) 


-0.02194 ± 0.00099 


0.03942 ± 0.00066 


0.7 


/2(1270) 


-0.699 ±0.018 


0.387 ±0.018 


2.1 


p(1450) 


0.253 ±0.038 


0.036 ± 0.055 


0.1 


Non-res 


-0.99 ±0.19 


3.82 ±0.13 


8.5 


/o(980) 


0.4465 ± 0.0057 


0.2572 ±0.0081 


6.4 


/o(1370) 


0.95 ±0.11 


-1.619 ±0.011 


2.0 


a 


1.28 ±0.02 


0.273 ± 0.024 


7.6 


a' 


0.290 ±0.010 


-0.0655 ± 0.0098 


0.9 



TABLE III; Complex amplitudes Urc'^'^^ and fit fractions of 
the different components {KsTr~, Ksi^^ , and tt+tt^ reso- 
nances) obtained from the fit of the D" -^ Kstv^tv^ Dalitz 
distribution from D*+ -^ D^tt^ events. Errors are statistical 
only. 



DALITZ PLOT ANALYSIS OF L»^ 



K+R- 



We study the decay Df — > K~^K 7r+ using a data 
sample of 240 fb~^. We focus particularly on the mea- 
surement of the relative decay rates „,„i ' ^ Z . ■. and 

B[Dt ^K'° (S92)K+) 



B{DJ^K+K--n + ) 



The decay Dj 



B{Dt^K+K-Tr+) 

0(lO2O)7r+ is fre- 



quently used as the D'^ reference decay mode. The im- 
provement in the measurements of these ratios is there- 
fore important. A previous Dalitz plot analysis of this 
decay used ^ 700 signal events [24]. We perform the 
present analysis using a number of signal events more 
than two orders of magnitude larger. 

We reconstruct the decay by fitting the three charged 
tracks in the event to a common vertex, requiring the x^ 
probability to be greater than 0.1%. We cleanly remove 
a small background from the decay £>*+ -^ D^+j^^n^ 
by requiring m^+K- < 1-85 GeV/c^. In Fig. [S] we show 
the invariant mass distribution of the reconstructed Df 



o 

o 

c 
o 
> 




0.5 0.75 1 1.25 

m(K"7:*) (GeV/c') 



0.75 1 1.25 

m(K";i*) (GeV/c^) 



0.75 1 1.25 

m(K"7:'") (GeV/c") 



1.25 
m(K";i*) (GeV/c^) 



FIG. 11: Legendre polynomials moments for the K^K (top) and K vr^ (bottom) channels of Df — > K^ K tt^. The dots 
with error bars are data points and the curves are derived from the fit functions. The results are preliminary. 



Decay Mode Decay fraction(%) 



Amplitude 



Phase (radians) 



K*{%%2yK+ 
0(lO2O)7r+ 

_/o (980)71+ 

ii'o*(1430)''A'+ 
/o(1710)7r+ 
/o(1370)7r+ 

K^{14ZQ)°K+ 
/2(1270)^+ 



48.7 ±0.2 ±1.6 
37.9 ±0.2 ±1.8 

35 ± 1 ± 14 

2.0 ±0.2 ±3.3 

2.0 ±0.1 ±1.0 

6.3 ±0.6 ±4.8 

0.17 ±0.05 ±0.30 

0.18 ±0.03 ±0.40 



1. (Fixed) 
1.081 ±0.006 ±0.049 

4.6 ±0.1 ±1.6 
1.07 ±0.06 ±0.73 
0.83 ±0.02 ±0.18 
1.74 ±0.09 ±1.05 
0.43 ± 0.05 ± 0.34 
0.40 ±0.04 ±0.35 



0. (Fixed) 

2.56 ±0.02 ±0.38 

-1.04 ±0.04 ±0.48 

-1.37 ±0.05 ±0.81 

-2.11 ±0.05 ±0.42 

-2.6 ±0.1 ±1.1 

-2.5 ±0.1 ±0.3 

0.3 ±0.2 ±0.5 



TABLE IV: The results obtained from the DJ — » K^K 7r+ Dalitz plot fit, listing fit-fractions, amplitudes and phases, 
errors are statistical and systematic, respectively. The results are preliminary. 



The 



candidate in the decay Df — > K^K^tt^ . For the Dalitz 
plot analysis, we use events in the ±2a mass window 
of the reconstructed Df candidate. We parametrize the 
incoherent background shape empirically using the events 
in the sidebands. In the signal region, we find 100850 
signal events with a purity of about 95%. 



The Dahtz plot for 
is shown in Fig. [TOl 



the _D+ -^ K^K^TT^ events 
m I'lg. IUJ| In the K+K~ threshold re- 
gion, a strong (/)(1020) signal can be observed, together 
with a rather broad structure indicating the presence of 
the /o(980) and ao(980) 5"- wave resonances. A strong 
ii'*''(890) signal can also be seen. We perform an un- 
binned maximum likelihood fit to determine the relative 
amplitudes and phases of intermediate resonant and non- 
resonant states. The complex amplitude coefficient for 
each of the contributing states is measured with respect 
to K*^K^. We summarize the fit results in Table HVl 
showing fit-fractions, amplitudes, and phases of the con- 



tributing resonances. The projections of the Dalitz plot 
variables in data and the ones from the fit results are 
shown in Fig. [TOl Further tests on the fit quality can be 
estimated using Y^ angular moments. These moments 
are shown for the K^K~ and K~'n'^ channels in Fig.fTT] 
The agreement between the data and fit is excellent. We 
find a rather large contribution from the /o(980)7r+, but 
with a large systematic uncertainty due primarily to a 
poor knowledge of the shape parameters of /o(980) and 
higher /o states. 



From the fit-fraction values reported in Tabic IIVI we 
make the following preliminary measurements: 



B(Dt^K+K-TT+) ~ 

g(D + ->J?'"'(892)i<'+) 
B(Dt^K+K-Tr+) 



0.379 ± 0.002 (stat) ± 0.018 (syst), 
= 0.487 ± 0.002 (stat) ± 0.016 (syst). 



CONCLUSIONS 



we have studied the amphtudes of the decays D'^ 



K-K+TT°, D^ 



r+TT^ Z?" -^ K'in+TT^, and D^, 



K+K'-K+. Using £)0 -^ K~K+ Tr° DaHtz plot anal- 
ysis, we measure the strong phase difference between 
the D^ and Z?" decays to K*{892)+K- and then am- 
plitude ratio, which will be useful in the measurement of 
the CKM phase 7. We observe contributions from the 
Kn and K^K^ scalar and vector amplitudes, and an- 
alyze their angular moments. We find no evidence for 
charged k, nor for higher spin states. We also perform a 
partial- wave analysis of the K'K"^ system in a limited 
mass range. We measure the magnitudes and phases of 
the components of the D° -^ 7r"'"7r~7r° decay amplitude, 
which we use in constraining the CKM phase 7 using 
B^ -^ D^^+^-^oK^. We measure the amplitudes of the 
neutral D- meson decays to the A'"7r^7r+ final state and 
use the results as input in the measurement of 7 using the 



decay B^ -^ D 



(*) 



_K^. Finally we parametrize the 



amplitudes of the Df -^ K^K 7r+ Dalitz plot and per- 
form precision measurements of the relative decay rates 

B{Dt-'K+K~-K+) B(Dt^K+K-'iT+) ' 
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